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Abstract 

This paper studies a slow-fast system whose principal characteristic is that the slow manifold is given by 
the critical set of the cusp catastrophe. Our analysis consists of two main parts: first, we recall a formal 
normal form suitable for systems as the one studied here; afterwards, taking advantage of this normal form, 
we investigate the transition near the cusp singularity by means of the blow up technique. Our contribution 
relies heavily in the usage of normal form theory, allowing us to refine previous results. 
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1. Introduction 


A slow-fast system (SFS) is a singularly perturbed ordinary dilferential equation of the form 

X = f{x,Z,£) 

£z = g{x,z,e), 

where x G K™, z € M" are local coordinates and where £ > 0 is a small parameter. The over-dot denotes the 
derivative with respect to the time parameter t. Throughout this text, we assume that the functions / and g 
are of class C°°. In applications (e.g [25]), z{t) represents states or measurable quantities of a process while 
x{t) stands for control parameters. The parameter e models the difference of the rates of change between 
the variables z and x. That is why systems like (1) are often used to model phenomena with two time scales. 
Observe that the smaller e is, the faster z evolves with respect to x. Therefore we refer to x (resp. z) as the 
slow (resp. fast) variable. The time parameter t is known as the slow time. For £ 7 ^ 0, we can define a new 
time parameter r by the relation t = er. With this time reparametrization (1) can be written as 

x'=ef{x,z,e) 

/ / \ 1^1 
2 =g{x,z,e), 

where now the prime denotes the derivative with respect to the rescaled time parameter r, which we call 
the fast time. Since we consider only autonomous systems, we often omit to indicate the time dependence 
of the variables. In the rest of this document, we prefer to work with slow-fast systems presented as (2). 
Observe that as long as £ yf 0 and / is not identically zero, systems (1) and (2) are equivalent. A first 
approach to understand the qualitative behavior of slow-fast systems is to study the limit £ —>■ 0. The slow 
equation (I) restricted to £ = 0 reads as 


x = f{x,z,0) 

0 = g{x,z,0). 

A system of the form (3) is called constrained differential equation (ODE) [14, 24]. On the other hand, in 
the limit £ —>■ 0 , a system given by ( 2 ) becomes 

x' = 0 

z'=g{x,z,0), 

which is called the layer equation. Associated to both systems, (3) and (4), the slow manifold S is defined 
by 


S' = {(x, z) e X M" I g{x, z, 0) = 0} , 


( 5 ) 


which serves as the phase space of the ODE (3) and as the set of equilibrium points of the layer equation (4). 
In the latter context, it is useful to recall the concept of Normally Hyperbolic Invariant Manifold (NHIM). 

Definition 1.1 (Normally Hyperbolic Invariant Manifold). Consider a slow-fast system given by a vector 
field of the form 


d d 

Ae = £/(x,Z,£)— -\-g{x,Z,£) — . 


( 6 ) 


The associated slow (invariant) manifold S = {g{x, z, 0) = 0} is said to be normally hyperbolic if each point 
of S is a hyperbolic equilibrium point of Xq . 
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Figure 1: A schematic representation of the persistence of a NHIM under the perturbation of the corresponding vector field. 
S denotes the slow manifold. Left-above: 5 is a set of hyperbolic equilibrium points of the layer equation. Left-below: S is 
the phase space of the constrained equation. Right: since 5 is a NHIM, it persists as an invariant manifold Se under small 
perturbations of the vector field. 


NHIMs are relevant in the context of the geometric study of slow-fast systems, see for example [10]. It 
is known that compact NHIMs persist under small perturbation of the vector field [15, 16]. In the 
particular context presented above, a normally hyperbolic compact subset of the slow manifold S persists as 
an invariant manifold of the slow-fast system Xg. We show in fig. I a schematic of the previous description. 
After this introduction, we turn into the subject of this paper. Our goal is to understand the dynamics of 
a particular slow-fast system which has one fast and two slow variables given as 

Ve =s(l+sfi)- -- (z^ + X2Z + Xi + sfs) —, (7) 

UXi 0X2 OZ 

where the functions fi = fi{xi,X 2 , z), for i = 1,2, 3, are smooth and vanish at the origin. The corresponding 
slow manifold is defined by 


S = {{xi,X2, z) e I -I- X2Z -I- xi = 0} . 


( 8 ) 


Remark 1.1. The slow manifold S can he regarded as the critical set of the cusp (or A^) catastrophe, which 
is given as [1, 5] 


V{xi,X2,z) 



i- 2 

-a:2Z -l-criz. 


(9) 


We denote by A the set of points in S at which S is tangent to the fast direction, that is 

X = [{X2, Z) & S \‘iz^ + X2 = H] ■ (10) 

In other words, A is the set of degenerate critical points of (9). See figure fig. 2 for a description of the slow 
manifold and the set A. 


S 



Figure 2: The manifold S is two dimensional and can be defined as the critical set of the potential function , X 2 , 2 :) = 

^ 2 :^ -|- ^X 2 Z^ H- xiz. The curve A is defined by the set of degenerate critical points of V. Geometrically, B is the set of point 
of S where S is tangent to the fast direction, and C denotes the cusp point. 

Our interest in studying (7) is due to the fact that the origin (a;i,a; 2 ,z) = (0,0,0) is a non-hyperbolic 
eguilibrium point of Xq. This implies that a compact subset, around the origin, of the slow manifold S is 
not a NHIM of Xq, and therefore, the Geometric Singular Perturbation Theory [10, 15, 16] is not enough. 
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1.1. Motivation 

There have been several studies, e.g. [18, 19], dealing with a SFS of the form 

Xg = £(1 + fi)~^ -+ xi + eh) —, (11) 

whose slow manifold is the critical set of the fold catastrophe. The next natural step is to consider the 
following case in the Thom list [22], i.e., a slow-fast system induced by the cusp catastrophe. That is 

Xe = £(1 + fi)^ -f -+ 2^2^ + xi + e/a) —. (12) 

OXI 0X2 OZ 

In [4], the system (12) is studied in a qualitative way. Here, however, we aim to refine the results by 
heavily using techniques from normal form theory. Moreover, we remark that the methods presented here 
are applicable to a larger class of slow-fast system given by 




e(l + /i) 


d 

dxi 


fc-i 




JL 

dxi 


k-l 

i=i 


rl-l 



which is called (regular) Hfe-SFS, see [13]. 

1.2. Statement 
We shall study the SFS 


(13) 


Xe = £(1 + /l)^-f -+ X2Z + Xi + e/ 3 ) —, (14) 

OXl 0 X 2 oz 

where the functions fi = fi{xi,X 2 , z,e) are smooth. To avoid working with an e-parameter family of vector 
fields as (14), it is customary to extend (14) by adding the trivial equation e' = 0, and thus consider a 
smooth vector field in which reads as 

00 fi fi 

X = e(l -I- fi)j^ -1- e/ 21 ^- {z^ + X 2 Z -I- xi -I- e/ 3 ) ——h 0 — . (15) 

0X1 0X2 OZ OS 

We regard (15) as a perturbation of “the principal part” F which is given as 

fi fi fi fi 

F = e- -hO^- {z^+ X 2 Z + X 1 ) —+0—. (16) 

uXi 0X2 OZ OS 

Note that in a qualitative sense, F contains the essential elements of X. To state our main result, we first 
define the sections 


E = [{xi,X2,z,e) &W^\xi = -x\] 
E“ = |(a;i, X2, z,e) G M"* | xi = a;{| , 


(17) 


where x\ > 0 and x{ > 0 are arbitrarily large constants. For £ > 0 but sufhciently small, the sections E“ 
and are transversal to the flow of X^. Next, let H : E“ —)■ E+ be the Poincare map induced by the flow 
of X^. We shall prove the following. 
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S- 



Figure 3: Description of our main result. We may choose appropriate coordinates at the sections and under which the 
invariant manifold Se is given by Z = 0. Moreover form (16) we have that all other trajectories starting at are exponentially 
attracted to the invariant manifold Se. In this paper we provide quantitative information regarding this exponential contraction. 


Transition along the cusp (see theorem 3.1). Consider a slow-fast system given by (15). Let S”, S+ 
and n : S” —S'*" be defined as above. Then, we can choose coordinates in and in S’*" such that the 
map n reads as 


niX2,Z,e) = iX2,Z,e), (18) 

where X 2 = X 2 + H{X 2 , e) (with H flat at {X 2 , e) = (0,0)), e = e and where 

Z ^ ^X 2 ,e) + Zexp (^-^{A{X 2 ,e) + e'l>{X 2 , Z,e))^ , (19) 

where A{X 2 , 0) > 0. Details of the functions A, and 'I' are given in theorem 3.1. In an heuristic way, this 
result is described in fig. 3. 

1.3. Idea of the proof 

Our proof consists of two main steps. 

1. From [12], it is known that there exists a formal transformation bringing (15) into 

F = e-^-h 0-^-( 2 ^ + X 2 Z -I- Xi) ——h 0—. (20) 

UXi UX2 UZ OS 

Then, by Borel’s lemma [5], the vector field F can be realized as a smooth normal form X^ = F + R 
of (15) and where R is flat at (a;i, 2 : 2 , z, e) = (0,0, 0,0). See more details in section 2.2. 

2. Based the previous normalization, next we use the geometric desingularization or blow up method (as 
introduced in [9]) to study the flow of the normal form X^ = F -\- R. This is detailed in section 3. 

Remark 1.2. With this document we aim at two goals: 

1. To refine the results of [4]. This is, we do not only provide a qualitative description of the transition 
n, but details on the differentiability of such a map is also presented. 
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2. To prepare a framework for the geometric desingularization of Ak slow-fast systems. These are a 
generalization of (15) given as 


X — e{l + /i) 


d 

dxi 


fc-i 




JL 

dx^ 


k-i \ 

z'^ +e/fc j 


— n— 
dz^^de' 


( 21 ) 


The rest of this document is arranged as follows: in section 2 we provide a brief recollection of preliminary 
results that will simplify our later studies. Next, in section 3 we pose our result and prove it by means 
of the geometric desingularization method and the results of section 2. For readability purposes, many 
technicalities have been put in the appendix. 


2. Preliminaries of slow-fast systems 

In this section, we provide a number preliminary results that will be used later in section 3. First of all, we 
consider slow-fast systems along normally hyperbolic regions of the slow manifold. Afterwards, we recall a 
result from [12] dealing with the normal form of (15). We remark that we only consider SFS with one fast 
variable. Let us be more precise with the type of SFS that we shall study first. 

Definition 2.1. A slow-fast system is said to be (locally) regular around a point po, if its corresponding 
slow manifold is normally hyperbolic in a some neighborhood ofpo- 

2.1. The slow vector field 

Let us consider a slow-fast system given by 

(9 d 

Ae = ^£/,(a:, 2 :,£)—-biL(a;,z,£) —, (22) 

i —1 * 

where x G z G M, and as usual 0 < e <C 1. Furthermore, assume that /(O, 0,0) ^ 0, 0, 0) = 0 and 

^(0, 0, 0) < 0. Thus Xg is regular around 0 € The slow manifold associated to (22) is defined by 

5'= {(a;,z) GM’”+i|iL(a;,z,0) =0}. (23) 

From the defining assumptions of (22), we have that S' is a NHIM in a neighborhood of the origin. By 
looking at the Jacobian of Ag at 0, it follows that there exists an m-l-1 dimensional a center manifold. Since 
X is smooth, we can choose a center manifold W for any £ < oo. The manifold W is given as a graph 
z = (j){x, e) where (p is function. 

Remark 2.1. Along the rest of the document we frequently make use of a finite class of differentiability. 
As it is customary in the present context, when we say that a manifold (or a map) is C^, we mean that such 
a manifold (or map) is £-differentiable for £ as large as necessary. 

The slow manifold S is naturally given by the restriction |£=o = S. Next, let us consider the vector field 
^Xg{x, (f>, e). Since W is locally invariant, it follows that ^X^ is tangent to W . Therefore the vector field 

A^'°“ = limiAg(cr,<^,£), (24) 

is tangent to S at each point of S, and we call it the .slow vector field. We remark that the slow vector field 

Xsiow jg defined whenever (p is invertible. 
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2.1.1. The slow divergence integral 

Associated to a regular slow-fast system and the corresponding slow vector field, the slow divergence integral 
is defined here. For this, let and be two sections which are transversal to the flow of given by 
(22). For e 7 ^ 0 but sufficiently small, these sections are also transversal to the slow manifold S. Let be 
a solution curve of chosen along a center manifold W , thus je is transversal to the sections 11 “ and 
In the limit e = 0, the curve 7 o is a curve along the slow manifold S. The idea now is to borrow 
the well-known divergence theorem [21] to get some sense on how the trajectories of X^ are attracted to S 
(recall that we made the assumption ^ < 0). The divergence of X^ (given by (22)) reads as 

divX,= +0{e). (25) 

oz 

We can now take the integral of divX^ along the orbit 7 ^ of parametrized by the fast time r, we have 

J divXg dr = J ^ o{e)^ dr. (26) 

The slow divergence integral is defined by 


m = 


div Xq dt^ 


'7o 


(27) 


where t is the slow time defined by the slow vector field Our goal then is to relate the divergence 

integral (26) with I. 

Proposition 2.1. Under the assumptions made in this section, we have that 


f div Xs dr =-{I{t) + o{l)), (28) 

where I{t) is the slow divergence integral. 

Proof. Recall that the slow vector field reads as = lim£_,.o -Xg{x,(j),s), where f = (j){x,e) is a 

function. By our assumptions, the curve is transversal to the sections and S"*" for e small enough. 
Without loss of generality we can assume that ye is parametrized by xi. Then let xf and x^ be defined by 
'leixf) = ye n and ye(a^7 = ye H E''". Next, the integral of the divergence of X^ along ye from E“ to 
E+ reads as 


div Xg dr = - 

e . 


dH{x, z, 0) 
dz 


+ 0{e) 


dxi 


fi{x,z,0) -ko(l) 


1 

£ 



dH{x,z,0) dxi 
dz fi{x,z,d) 



1 

e 


div ATo dt + o(l) 


where t is the slow time induced by which in coordinates means that 7 ^ = /i. 


(29) 


□ 


Observe that the slow divergence integral is a first order approximation of the divergence along orbits of Xg. 
This will be useful when presenting our main result in section 3. 


7 



2.1.2. Normal form and transition of a regular slow-fast system 

Now we consider the problem of finding a suitable normal form of a regular SFS. The following is a well- 
known result but we recall it here for completeness. 

Proposition 2 . 2 . Consider a regular slow-fast system on given by 

i=i ■' 

where {u,vi,... ,Vm, z,e) € where the functions f\ = fi{u,v, z,e) and gj = gj{u,v, z,e), for 2 > 

j > k — 1, are smooth and where the function H = H{u,v, z,e) is smooth with i/(0,0,0,0) = 0 and 
^(0,0, 0,0) < 0. Then, the vector field X is -equivalent to a normal form given by 


X, 


N 


d d ^ d 


i=i 


dV, 




(31) 


where {Z = 0} corresponds to a choice of the center manifold W of Xg. 

Proof of proposition 2.2. The first step is to divide the vector field X hy 1 -\- fi. In a sufficiently small 
neighborhood of the origin this is a smooth equivalence relation. That is F = reads as 


^ d ^ ^ d ~d 


i=i 


dvn 


dz 


(32) 


where gj, for 2 > j > A: — 1, and T1 are smooth with H{0) = 0 and ^(0) < 0. Now we note that the 
origin of is a semyhyperbolic equilibrium point with {u, v, s) being center coordinates and 2 being the 

hyperbolic coordinate. We can now use Takens-Bonckaert results on normal forms of partially hyperbolic 
vector fields [2, 3, 23]. Thus, there exists a change of coordinates (maybe respecting some constraints if 
required) under which Y is conjugated to 


Tfl »-\ 


i=i 


(33) 


where Gj = Gj{U,V,s), for 2 > j > k — 1, and H = H{U,V,e) are functions, and where {Z = 0} 
corresponds to a choice center manifold which we denote by W . We remark that in the vector field Y, the 
functions Gj and H are independent of Z. Furthermore we have 

H{0, 0,0) = — (0,0,0,0) < 0. (34) 

This means that in a small neighborhood of the origin Y can be divided by \H\. In other words, Y is 
C^-equivalent to 


r\ Tn r\ r\ 




'dY 


dZ' 


(35) 


where 11(0, 0, 0) 7 ^ 0 and Kj = Kj{U, F,e), for 2 > j > A: — I, are . Next, since W = {Z = 0} is invariant 
under the flow of y, we can study the restriction y\z=o. This is 


p. m p. 

y\z=, = eQ^,+Y.eK,—. 


dU 


(36) 



For e 7 ^ 0, the vector field 3^|z=o is regular because tj(0, 0, 0) ^ 0. Thus, by the flow-box theorem, there 
exists a change of coordinates, depending in a way on e, under which ylz=o can be written as 


d d 


This implies that y is C^-equivalent to 


r\ r\ 


i=i 


dVi dZ' 


as stated in the proposition. 


Motivated by proposition 2.2 let us now discuss the dynamics of the vector field 

r\ 

__I Vn—-7— 


i=i 


dVi dz' 


The slow manifold S, corresponding to the normal form (39), is given by 

S = {e = Q, Z = 0} . 


(37) 


(38) 


□ 


(39) 


(40) 


Furthermore, we can parametrize the solution of (39) by U. Let us define the sections 

j:- = {{U,V,Z,e) eRxW^ xRxR\U = U-} 
s+ = {([/,F,Z,£) e M X M™ X M X M|{7 = [/+} , 


where U < C/+. The sections E and E+ are transversal to the manifold S and therefore, for e 7 ^= 0, are 
also transversal to the flow of (39). Associated to these sections, we define the transition 

n : E" ^ E+ 

. , (42) 

{V,Z,e)^iV,Z,e). 


To compute the component Z we only need to integrate = — ^Z. Then it follows that Z = Z{T), where 
T is the time to go from E“ to E+, which is T = Uf — Ui. Then it follows that 


V = V 

Z = Zexp(^-^{Uf-Ui)^ (43) 

£ = £. 

Observe the particular format of the transition If. The Z component is an exponential contraction towards 
the center manifold {Z = 0}. Maps with this characteristic appear frequently in our text and also in several 
other cases where slow-fast systems are studied. Therefore, in appendix A we discuss in a rather general 
way, the properties of such maps. 

2.2. Formal normal form of Ak slow-fast systems 

In this section we recall a normal form of the so-called Ak slow-fast systems. A proof can be found in [12]. 
This normalization is important since it eliminates many unwanted terms from the system being studied 
here. 
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Definition 2.2. Let /c G N with k >2. An slow-fast system (Ak-SFS) is an ODE of the form 

x[ = e(l + /i) 
x' = efj 

k-1 

e' = 0, 

where j = 2,... ,k — 1, and where the functions fi = fi{xi, ..., Xk-i, z, e), for 1 < i < k, are smooth. 

Remark 2.2. 

• The system investigated in this work is an A^-SFS. 

• The slow manifold associated to an A^-SES is defined by 

k-1 'I 

(x, z) G I ^ = 0 I . (45) 

i=l J 

The manifold S can equivalently he defined as the critical set of an Ak catastrophe [1 ]. Hence the name 
Ak-SFS. 

Locally, we can regard (44) a.s X = F -{- P where F and P are smooth vector fields of the form 





and 



k-1 




d 

dxj 


d 


(46) 


k-1 




d ^ . d 

Z, - 1 - ^fk-^ 

OXi OZ 


de 


(47) 


respectively and where g = — \ z^ + Y^\=i We refer to F as the “principal part” and to P as 

the “perturbation”. Briefly speaking we want to eliminate, via a change of coordinates, the perturbation. 
The procedure of normalizing the vector field X is motivated by [20], where normal forms of analytic 
perturbations of quasihomogeneous vector fields are investigated. The relevant result is the following 

Theorem 2.1 (Formal normal form [12]). Let k>2 and let X = F -\- P be a smooth vector field where 


and where 


F = e 


dx\ 


k — l a 

* 


- 2 : 


k-1 

i=i 


•i-i 



k-1 


p = Y.p^ 

i=l 


d_ 

dxi 


+ Pk 


d_ 

dz 


+ 0 


de' 


(48) 


(49) 


where each Pi = Pi{xi, ..., Xfe_i, z, e) is a smooth function. Assume that the following conditions are satisfied 


1. Pj(a;i,...,a;fc_i,z,0) = 0, 
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p{Pi) > 2fc — i + 1, 

where Pi denotes the Taylor expansion of Pi and p{Pi) is the quasihomogeneous order of the polynomial Pi. 
Then, there exists a formal diffeomorphism such that = F- 

In words, theorem 2.1 shows that X and F are conjugated via $. It follows that, by Borel’s lemma [5], 
the formal vector field = F can be realized as a smooth vector field X^ = F + P where P is flat 
at {x,z,e) = (0,0,0). This has important consequences in the geometric desingularization of an A 3 -SFS, 
presented in the following section. 


3. Geometric desingularization of a slo-w-fast system near a cusp singularity 


In this section we study an A 3 slow-fast system based on: a) the techniques introduced in section 2 and in 
appendix A, and b) the blow up method. To simplify the notation, let us now write the A 3 -SFS as 

^ + a + e/ 3 )^ + 0^, (50) 

where thanks to theorem 2.1, the smooth functions fl = fi{a, b, z, e) are flat at the origin of We invetigate 
the transition associated to (50) between the sections 


S = {(a, 6, 2 , e) £ I a = —a , 0 > 0} 

= {(a, 6, 2 , e) £ I a = a'*', z < O} , 


(51) 


where a >0 and o’*" > 0 are arbitrarily large constants. However, since the trajectories of X spend a long 
time along regular parts of S, it will be useful to define the “entry” and “exit” sections 


S®" = {(a, b,z,e) £ | a = —oq, z > O} 

= I (a, 5, z, e) £ K'* | a = oq, z < 0} , 


where oq is a positive but sufficiently small constant, for reference see fig. 4. 



Figure 4: Qualitative representation of the investigation performed in this section. The sections and are arbitrarily 
close to the cusp point. On the other hand the sections and S'*" (not shown) are parallel to and but far away from 
the cusp point. In a qualitative sense, we will construct an invariant manifold Adg and then extend it all the the way up to 
the sections S“ and S"^. Our analysis aims for simplicity and thus depends extensively on the usage of normal forms. This, 
of course, makes our results coordinate-dependant. 


It will be clear from our analysis in the blow up space section 3.2 that the section E needs to be partitioned 
as follows. 
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Definition 3.1 (The inner layer and the lateral regions). Let 0<L<M<oobe constants. The inner 
layer C is defined as 

E- 3 E“‘“'' = {(6,z,e) e E- 116 | < . (53) 

On the other hand, the lateral regions are defined as 

E- D E+^ = I(6,0, e) G E- I 6 > Le 2 / 5 | 

. ( 54 ) 

E- D E-^ = |(6,z,£:) G E- I - 6 > 

Note that the set |E“"®'', E"*'*’, E”**} is an open cover ofY,~, see fig. 5 


pinner 



Figure 5: The section E“ needs to be partitioned into three subsections: the inner layer and the lateral regions 

From a qualitative point of view, these three layers correspond to three different types of trajectories: 1. Trajectories 
starting at pass close to the cusp point. Observe that = {6 = 0} and then corresponds to a solution of 

the associated CDE passing exactly through the cusp point. 2. Trajectories starting at pass sufficiently away from the 
cusp point along the regular side of the manifold S. 3. Trajectories starting at pass sufficiently away from the cusp point 
along the folded side of the manifold S. 


We are now in position to present our main result. In the following theorem, we characterize the transition 
n : E“ —>■ E+ under a suitable choice of coordinates at the section E“ and E+. Furthermore, we give details 
on the differentiability of this map according to the cover of E“, see definition 3.1. 

Theorem 3.1 (Transition map of an A 3 -SFS). Let X he an A 3 slow-fast system. This is, X is a vector 
field defined by 

^ = + + e/ 2 ^ - + 6 ^ + a + 6 / 3 ) ^ + 0^, (55) 

where each fi = fi(a,b, z,e), i = 1,2,3, is smooth. Let the sections E“, E+ he defined as above. Then we 
can choose suitable -coordinates {B, Z, e) in E“ and -coordinates {B, Z, e) in E’*' such that the transition 
n : {B, Z, e) I—)■ {B, Z, e) is an exponential type map of the form 


n{B,z,s) 


^B + h, (j>{B, e) Z exp 


A{B,e) + -^{B,Z,e) ^ 


(56) 


where h is flat at the origin, A > 0 is , (j) is -admissible with (f>{B,0) = 0, and d' is C^-admissible 
with '^{B,Z,0) = 0, see appendix A for the definition of -admissible. Moreover, we have the following 
properties of the function A, (j) and . 
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1. —A{B,0) = 1(B) where I is the slow divergence integral associated to (55). 

2. Restricted to (B,Z,e) G there are functions 4> and such that 


'it(B,Z,e) = Zj , 

where 4> and ^ are -functions with respect to monomials (see definition A.2) with p = . Note 

that in this domain, p is well defined in the sense that p is bounded by a constant os £ —>■ 0 . 

3. Restricted to (B,Z,e) G 1]+^, there is a function ^ such that 

(/)( B ,£)=0 

'i>iB,Z,e) = § (|i?|l/^£l/^£ln(|i?|),a,z) , 


where ^ is a -function with respect to monomials (see definition A.2) with a = e\B\ Note that 
in this domain, a is well defined since \B\ > 0. 

4- Restricted to (B,Z,e) G 1]“^, there are functions (p and ^ such that 


(p(B,e) = f (\B\^/^,a^ 

^iB,Z,e) = ^{\B\^/\e^/^e\n(\B\),a), 


where (j) and ^ are -functions with respect to monomials (see definition A.2) with a = s\B\ 
Note that in this domain, a is well defined since \B\ > 0. 


Sketch of the proof. The first step is to recall theorem 2.1, which shows that X is formally conjugate to 

F = e^+0^-(z^ + bz + a)^+0^. (60) 

oa ob ' ' oz oe 

Next, by means of the Borel’s lemma [5], the vector field F can be realized as a smooth vector field 
= F eFI where H is flat at (a, b, z, s) = (0,0, 0, 0). Thus, from now on, we only treat an ^ 3 -SFS given 
as 


X — £(1 + ^/i)^ + ~ A bz + a + sf^^ — + 0^, (61) 

where each fi = fi{a, b, z, e) is flat at (o, b, z, e) = ( 0 , 0 , 0 , 0 ). 

Another important ingredient of the proof is the blow up technique, which is described in section 3.1. This 
method provides several local vector fields whose corresponding transitions are of exponential type, refer 
to appendix A. Later all these local transitions are composed to produce an exponential type transition 
between the sections and E"*". Along the analysis of the local vector fields (in the blow up space) we 
will take advantage of the flatness of the higher order terms of X. The complete proof follows sections 3.1 
to 3.5 and is given in section 3.6. 

Now, assuming that the transition 11 is of the form (56), we can show that A(B,0) is given by the slow 
divergence integral of X. For this, let us recall the Poincare-Leontovich-Sotomayor formula [7], which in 
general is given as follows. 
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Proposition 3.1. Let X be a veetor field on a manifold M" with a volume form LI. Let E and E'*' be two 
open seetions of M and transverse to the flow of X. Let 7 ^ be an orbit of X along a center manifold W 
of X, starting at p = js L E“ and reaching g = 7 e n E+ in finite time. Let 11 : E“ —)■ E+ he the transition 
map defined in a neighborhood of p. If ip~ : U —)■ E“ and : V —)■ E+, with U C and V C 

are coordinates in E“ and in E+ respectively, then 


det (£> ((V'"'’) ^ollo'i/) )) (s ) 


{n{p),Df {s )xX{p)) 
{n{q),Df+is+)xX{p)) 


divn X dr 


(62) 


where s = {ip ) ^{p) and s^ = {ip'^) ^{q). The integral is taken along the orbit from p to q parametrized 
by the fast time t. 


So we have the following. 

Proposition 3.2. Consider an A 3 -SFS and assume that the transition If : E“ —>■ E+ is given by (56). 
Then —A{B,0) = I{B), where I{B) is the slow divergence integral associated to the A^-SFS. 

Proof. The only relevant component is Z, so denote by II^ the Z-component of 11. The factor multiplying 
the exponential in (62) can be taken as a constant C > 0. Then we have that (62) for the vector field of 
theorem 3.1 reads as 


dIVz 

dZ 



divo X dr ] . 


(63) 


Using the properties of the slow divergence integral described in section 2.1.1, and since C 7 ^ 0, we have 


dAz 

dZ 


= C exp 


divQ X dr 


7e 


exp div Jfo dt + elnC + 0 ( 1 )^ ^ 


(64) 


'70 

= exp ( ^ (/ + 0 {e)) ) , 


where / is the slow divergence integral of X along a curve in the slow manifold S from E“ to E'*'. In 
principle, the limit e —> 0 of (64) is not well defined. However, according to our theorem 3.1, we have by 
differentiating (56) w.r.t. Z 


dliz f A{B,e)+e'i>{B,Z,e)\ 

-dz=^^^\ -^- )■ 

Identifying (64) with (65) and taking the limit e —>■ 0 we have indeed that 

lim(/ + 0{e)) = lim(—H(H,e) + e'I'(H, Z,e)), 

£—>■0 £—^0 

which shows the claim. Note that the slow divergence integral in the coordinates (a, b, z) reads as 

/(6) = j(6,c+)-/(6,r), 

where straightforward computations show that 

/>,C) = ^C" + 2C"fc + 6"C, 

and where is a constant defined by (a^, 5, C^) € E=*= n S. 


(65) 


( 66 ) 


(67) 


( 68 ) 
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On the other hand, in normal coordinates and along regular parts of the slow manifold, the ^ 3 -SFS can be 
written as (see section 2 . 1 . 2 ) 


X(.4.B,Z..)=4+oA-^A+o|:. (69) 

In these coordinates the slow divergence integral reads as 

I = A+-A-, (70) 

where A'^ and A~ are the corresponding parametrizations of and (respectively) in the coordinates 
(A, B, Z, e). 

□ 

3.1. Blow-up and charts 

Let us briefly recall the blow up technique, for more details see e.g. [ 8 , 9, 17]. The vector field X (50) is 
quasihomogeneous [1, 12]. Therefore, it is convenient to use the quasihomogeneous blow up. This technique 
consists on performing a coordinate transformation defined by 

a = r^d, b = r'^b, z = rz, e = r^e, (71) 

which is called the blow up map, and where = 1 and r € [0, + 00 ). That is (a, b, z, e, r) € 

X IR+. Since e > 0, we can restrict the coordinates to £ > 0. Note that x {0} is mapped, via the blow 
up map (71), to the origin of The powers or weights of the blow up map (71) are obtained from the 
type of quasihomogeneity of X. 

Let us denote by $(a, b, z, e) the blow up map (71). This map induces a smooth vector held X on x IR+ 
dehned by = X. It is often the case in which the vector held X is degenerate along x {0}. Then one 
dehnes another vector held X by X = 4 X for a well chosen positive integer m so that X is non-degenerate 
along X {0}. Since r € M"*", the phase portraits of X and X are equivalent outside x {0}, and therefore 
it is equally useful to study X instead of X. One obtains a complete description of the local how of X near 
the the cusp point by studying the how of X for (a, 6, z, £, r) G x [0, rg) with rg > 0 sufhciently small. 
For problems of dimension greater than 2, performing computations in spherical coordinates becomes tedious. 
Therefore, it is more convenient to consider charts which parametrize hemispheres of the ball x [0,rg). 
In the present context, the useful charts are 

Xen = {a = -1} , = {o = 1} , K, = {s = 1} , K± = {b = ±l} (72) 

and we always keep r € [0,rg). The previous setting is also known as directional blow up. A qualitative 
picture of the charts is given in hg. 6 . 

Briehy speaking, our analysis goes as follows: hrst, we perform a local analysis on each chart given in (72). 
Next, we compose (“glue”) the local results to provide a full description of the how of X (50) in a small 
neighborhood of the cusp point. In this way, we construct an invariant manifold from E®" to E®’^. Later we 
“push away” this invariant manifold all the way up to the sections E“ and E’*' along regular parts of the 
slow manifold S. 

To avoid confusion of the coordinates we adopt the following notation. Any object O dehned in the chart 
Ken is denoted by Oi. Similarly any object dehned in the chart Kex is denoted by O 3 . Finally, an object O 
dehned in either of the charts K^ or K± is denoted by 02 - 

3.2. Analysis in the chart Ken 

Taking into account our notation convention, the blow-up map in this chart is given by 

a = -rf, b = rjbi, z = rfzi, e = rlei. (73) 
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Figure 6: The blow up space and the charts. Each chart Ki parametrizes a region of the ball x [0, ro). A local analysis in 
the charts provides a full picture of the dynamics of the vector field X. 


The corresponding vector field in this chart (after multiplication by 3) has the form 


r[ = -eiri 


b[ = 2£ibi 


(i+/.) 




z[ =-3 (z? + 6 izi - 1 - i£i 2 i) + r^ei/s 

;'i = dsf (l + fij 


where the functions fi = fi(ri, bi, zi, Si) are flat along n = 0 , recall that x {r = 0 } 0 G via the 

blow up map. We study a transition IIi : A®" —>■ Ay where 

= {(ri,bi,zu£i) G In = ro,£i < S, zi > O} 

r ^ ^ (70) 

= {(n,bi,zi,£i) G I £i = (5,n < To}, 
where ro and S are sufficiently small positive constants. 

Remark 3.1. The section Ay corresponds to S®” in the blow-up space, that is S®" = $(Ay), where $ 
is the blow-up map (73). This implies that trajectories of X crossing S®" correspond to trajectories of 
crossing Ay. 

Before going any further, let us provide a qualitative description of Agn as in [4]. This process can be 
repeated, following similar arguments, in all the local charts; however, for brevity we only detail it for the 
current one. 

Qualitative description of the flow of Agn. The subspaces {ri = 0}, {ei = 0} and {ri = 0} n {si = 0} are 
invariant. Therefore, it is useful to study the flow of Ag^ restricted to the aforementioned subspaces. 

Restriction to {ri = 0} n {si = 0}. In this space Agn is reduced to 

b[ =0 


The set 


z} = -3 {zf + bizi - l) . 


7 i = {{bi,zi) I zf + biZi -1 = 0 } 


is a curve of equilibrium points. The phase portrait of (76) is shown in figure 7. 
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Figure 7: The phase portrait of Xen restricted to the invariant space {ri = Olfllei = 0}. The shown curve is 71 and it comprises 
a set of equilibrium points. Note that locally, all trajectories with initial condition 21 (0) > 0 are attracted to 7 i|{ 2 ^> 0 }- 


Remark 3.2. All the trajectories of (76) restricted to an initial condition Zq > 0 are attracted to the curve 
7 i|zi>o- Furthermore, due to our definition of Af^, we are interested only in trajectories satisfying this 
initial condition. Thus, from now on, we restrict our analysis to the subspace {zi > 0}. 

Restriction to {ei = 0}. In this space is reduced to 

r[ = 0 

b[ = 0 (78) 

7 = -3 (7 + bizi - 1) . 

The set Ti = {(ri, 6 i , zi) I zf + &i 2 i — 1 = O} is a surface of equilibrium points given by Ti = (ri, 71 ). Since 
r[ = 0, the phase space of (78) is foliated by two dimensional leaves in which the flow looks like fig. 7. 

Restriction to {ri =0}. In this space is reduced to 

A 

e'l 

Once again, the set 71 = {( 6 i,zi,£i) | £1 = 0, z > 0, zf + &iZi — 1 = 0} is a curve of equilibrium points. 
The Jacobian of (79) evaluated along 71 shows that, for small enough £ 1 , there exists an invariant center 
manifold that passes through 71 . Furthermore, the non-zero eigenvalue corresponding to the z-direction is 
negative along 71 . The phase portrait of (79) is shown in figure 8 . 

Observe that the bi and the £1 directions are expanding. It is important to know the relation between such 
two expanding variables. We have 


— 2£i6i 


1 


= -3 zf -b 61Z1 - 1 - -£iZi 


(79) 


5e?, 


db\ 2 h\ 
dei 5 £1 ’ 


which has the solution 


hi = h\ 



(80) 


(81) 
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Figure 8: Phase portrait of (79) restricted to zi > 0. The shown surface is an invariant center manifold, which is attracting in 
the 2 :i-direction. 


where < bi and e* < ei are the initial conditions, that is (6J,£j) = (&i,£i)|a 5 "- It is important to look 
at the ratio of initial conditions \/5 ■ This ratio tells us that bi is bounded as ei —>■ 0 (and therefore as 


el —>■ 0) if and only if G O ((ep^'^^). In other words, if the initial condition is not of order 0((ep^/®) 
then the value of bi at A®’' blows up as e* —>■ 0. This leads us to partition the section A™ into three open 
regions as follows. 


y^en,inner 
^en.bi 


Afl 

An 

An 


IbiKMeJ''® 


(82) 


where 0 < K < M < oo. Observe that the open sets ^en,bi 

A®". Accordingly, these sets induce an open cover of the entry section 
fig. 9 for a representation of the aforementioned partition. 


and A™’ form an open cover of 
E®'' via the blow up map (73). See 



Figure 9: Partition of A|“. Trajectories crossing through corresponding to the inner wedge area, have a continuation 

on the chart K^. On the other hand, outside we must consider the lateral regions and 

Based on the partition of the entry section A™, we define three transitions as follows 


j-^inner . ^en,inner _^ y^ex 

n+bi . ^en.+bi ^ ^ex,+bi (gg) 

n-br . ^en.-bi ^ ^ex,-bi ^ 


where 
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( 84 ) 


= {iri,bi,Zi,Si) e |ei = ^,ri < tq} , 

^ex,±6i ^ ^(^ri,bi,zi,si) e I 6i = ±77,ri < ro} . 

To finish with the qualitative description, note that there exists a (non-unique) 3-dimensional center manifold 
Wi , which is shown to exist by evaluating the Jacobian of X^n all along the surface 


Ti = {(ri,6i,zi,ei) | £1 = 0, 01 > 0 zf + bizi - 1 = O} . (85) 

Moreover, by the analysis provided above, the center manifold Wi\zi>o is attracting for ei small enough. 
Note that Wi |ei=o = hi. This means that Wi can be interpreted as a perturbation of the slow manifold S, 
written in the coordinates of the current chart. See fig. 10 for a representation of the previous exposition. 


£1 



Figure 10: Phase portrait of the trajectories of Xen depending on their initial condition. If the trajectories satisfy the estimate 
y E then they arrive to in finite time. If the estimate y E is not satisfied, then we must choose one of 

the outgoing sections in order to have a well defined transition map. 

Let US recall that the vector field Xen is of the form 

r[ = -eiri (^1 -h /i) 

X„: (l + A) + rfef/, ^ 

z[ =-S{zf + biZi-l-^eiZi)+rleif3 

e[ =5sl(l + fi'^ 

We now proceed to describe the transitions Hi given by (83). For this, first we write (86) in a suitable 
normal form. Next, based on this normal form, we compute the corresponding transition. 

First of all, let us move the origin to the point (ri,b,zi,ei) = (0,0,1,0). This is done by defining a new 
variable Ci by — 1. With this variable we have a new local vector field Fen which is defined by 
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( 87 ) 


r'l =-Siri(l + fi^ 

b[ = 2eibi (l + /i) + rlelf2 

e'l =5el{l + h) 

,Ci = ~3G(6i, £1, Cl) + £1^, 

where G(0,0,0) = 0 and f^(0,0, 0) = 3. Now, we want to write ICn in a suitable normal form. From 
proposition C.l, we know that ICn is equivalent to 




N 


{ r'l = -SiTi 
= 2 £iSi 
£' l = 5£? 

Z[ =-9{l + Hi{n,Bi,ei))Z 


1 ) 


( 88 ) 


where Hi is a G^-function vanishing at the origin. This normal form is convenient since the chosen 
center manifold Wi is now simply given by Wi = {Zi = 0}. Furthermore, from the format of X^, it is 
evident the “hyperbolic nature” of the flow restricted to the center manifold: the restriction of X^ to the 
center manifold Wi has a simple structure, namely 


{ r'l = -eiri 

B'l = 2siBi (89) 

£'i = 5el 

Note that for £i ^ 0, the vector field is hyperbolic. 

The vector field X^ is of the form studied in proposition C.4, therefore we have that the transition 

n““’' ■.{Bi,ei,zi)^ {ri,Bi,Zi) (90) 


is of the form 


f 61^/5 
n = TO [j) 


Bi=Bi{ - 

£l 


2/5 


Zi — Zi exp [ — -—(1 + aiGi \nsi H- SiGi) j , 


(91) 


where ai = Q;i(ro|i?i|^/^,ro£/^) and Gi = Gi(ro|i3i|^/^,ro£/^,^) where fi = BiSi Recall that for this 
transition we have the condition Bi € 0(£/^) so ^ is well defined. 


On the other hand, the transition 

: (i3i,£i,Zi) I— >■ (fi,£i,.Zi) 


(92) 
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is (see proposition C.4) of the form 


ri=rQ\ — 

V 


£l — £i 


Jl_ 

Bi 


1/2 


5/2 


Zi = Ziexp ( - —(1 +/3iei ln(|Bi|) + SiHi) ) , 


( 93 ) 


where /?i = /?i(ro|-Bi|^/^,roej^^) and Hi = _ffi(ro|i3i|^/^,ro£i^^,cr), where cr = ei|-Bi| Note that since 
Bi ^ 0(£^^^), a is well defined. We observe that the transitions and are exponential type maps. 


3.3. Analysis in the chart Kg 

Taking into account our notation convention, the blow-up map in this chart is given by 


a = r|a2, b = r|62, -z = rlz2, £ = rf. 


(94) 


Then, the blown up vector field reads as 


u =0 
I 02 =l + gi 

" ■ j 6 ' = r^g2 

U2 = - (^2 + ^2^2 + 02) + 53, 


(95) 


where the function gi = gi{r 2 , 02 , & 2 , Z 2 ) are flat along r 2 = 0. Note that in this chart r 2 acts as a parameter 
and that the flow is regular. Furthermore, note that Xg is not a slow-fast system, but a regular vector field. 
From the equation o^ = 1 -I- 51 , we define the following “entry” and “exit” sections. 

= {(^ 2 , 02 ,^ 2 ,2:2) I 02 = -^0, 22 > 0 } , 

(96) 

= {(^2, 02,^2,2:2) I 02 = 4 lo, Z2 < 0 }. 


Therefore, we define a transition HI as 


nf :A 


en,£ 

2 


^ A 


en,£ 

2 


{r2,b2,Z2) {f2,b2,h)- 


(97) 


Since (95) is regular, by the flow box theorem all trajectories starting at A™’^ arrive at in finite time. 
Moreover, the transition HI is a diffeomorphism and then, from (95) we have that HI reads as 


n2£(r2, b2,Z2) = if 2, b2,Z2) 

= {r2,b2 + hh 2 ,(/)i(r 2 , 62 ) + <(' 2 (^ 2 , &2)(1 + (I>3{'r2,b2, Z2))z2), 

where the <j)iS are smooth functions. Observe that in this chart, the transition is not an exponential type 
map. 


21 



3.4- Analysis in the chart i^ex 

Taking into account our notation convention, the blow-up map in this chart is given by 


a = rf, b = rlh, z = r^za, e = rfeg. 


( 99 ) 


Then, the blown up vector held reads as 


ra = eara (l + /i) 

=-2"3fc3(l + /i)+rfe^/2 ^ 

z' = -3 (^1 -h baZa + 1 + \eaza) + rlsafa 
4 =-5£i(l + /i) 

where the function fi = fi{ra,ba,S 3 , Za) are flat along ra = 0. Observe that the vector held Xex resembles 
the vector held Therefore, we have a similar behavior of the trajectories, the main difference is that 

in the case of Xex, there is one expanding (ra) and three contracting {ba, £3 and za) directions. The how of 
Xex is obtained following similar arguments as for the how of Xen- 

From the fact that Xg^ has three contracting and one expanding direction, we dehne the entry sections 

= {{i"3,ba,£3,za) ■■ £3=S, za< 0, ra < ro} 

= {(£3,63, £3,23) : ba = r ], Za < 0, ra < tq} (101) 

= {{ ra , ba ,£ a , za ) : 63 = -??, za < 0, rs < ro} , 
where all the constants are positive and sufficiently small, and the exit section 


Ar = {{r3,ba,£3,za) : ra = ro, za < 0, £3 < S, | 53 | < r]}. 
Then, accordingly, we dehne three transition maps as follows 


( 102 ) 


n^3 . ^en.e ^ ^ex 


: ira,ba,za) ^ ( 63 , £ 3 , 23 ) 


n+b3 . ^ ^ex 

: (£ 3 , £ 3 ,^ 3 ) (^3, £ 3 ,%) 

Ua^^ : ^ A®’' 

: (£3, £3,^3) ( b 3 ,£ 3 , Z 3 )- 


(103) 


Now we proceed to write Aex in a normal form just as we did with Xg^ in section 3.2. Following proposi¬ 
tion C.l we have that Agx is equivalent to 


{ £3 = £3£3 

B'a = -2eaBa 

£' = -5£i 

Z' =-9{l + Ha)Za, 


( 104 ) 
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where H 3 = 8 ^, 63 ) is a function vanishing at the origin. Just as in the chart Ken, there exists 

a three dimensional center manifold Wg associated to and which has been chosen such that Wg = 
{Zg = 0}. Since is the only expanding direction, we take as transition time Tg = In . This transition 

time is computed from the dynamics restricted to that is, from the equation rg = r In contrast to 
what happened in the chart Ken, the time T 3 is well defined for all the three transitions Ilg^, and 

Ilg Following proposition C.4 we have 


S 3 

£3 

^3 



l + a3e3 Inr3+£3K3 


(105) 


where = Q; 3 (r 3 |i? 3 |^/^,r 3 e 3 ^^) and J /3 = K 3 (r 3 |i 33 |^/^,r 3 ). Therefore, by taking the definitions of 
the entry sections we have 


n?(r3,B„Z3) = (b3 (5)'. i (())'. Z3exp (-4 { 0 ‘ - 

nf(x3,X3.Z3)=(±.3(;^) ,3,(5) .Z33 Xp(^-A(^(5) 

Observe that these transitions are of exponential type. 

3.5. Analysis in the charts K^i 
In this section we study the local flow at the charts and K_i. In a qualitative sense, these charts come 

into play when the initial condition 69 = does not satisfy the estimate 69 € This implies that 

the corresponding trajectory passes away from the cusp point. The chart K_|_j “sees” trajectories with initial 
condition 6 |Een > 0 while K_i “sees” trajectories with initial condition 5|sen < 0 . 


1 + a^5 In rs + JKg J J J 
- 1 + a3S3 lnr 3 + S3II3 


(106) 


Analysis in the chart K_^i 

In this chart the blow-up maps reads 


a = r\a 2 , b = r^, z = r 2 Z 2 , e = rle 2 . 
Then we have that the blow-up vector field is given by 


K 


+ b ■ 


2 ~ £ 2 /r 

4 = £ 2(1 + /aa) + £2^02 

2 = -£ 2 /e 2 

'2 = ~(^2 5 - Z2 + 02) + £2/2; 


where all the functions fe are flat along {r 2 = 0}. Observe that the set 


(107) 


(108) 


r2 — {(^2,02, £2, Z2) \ S2 — 0, Z2 + Z2 + a2 — 0} 


(109) 
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is a NHIM of However, X_^_l is not exactly a slow-fast system since £2 ¥" Oj but the restriction of X_^_l 
to {r 2 = 0} is indeed a slow-fast system. This restriction reads as 


^-|-bl{'’2=0} 


{ 02 = £2 

e'2 =0 

z'2 = -(^2 +-2^2 + 02). 


( 110 ) 


Remark 3.3. The subspace {r 2 = 0} is invariant. Moreover, since X_^_i is a flat perturbation o/X_|_j|{,r 2 =o}> 
it is equally useful to study the restriction X_|_^|{r 2 =o} • After all, by regular perturbation theory, their flows 
are equivalent. 

The slow manifold of .Y_|_h|{r- 2 =o} is defined by r 2 |r- 2 =o and is normally hyperbolic. Let us define the sections 


= {(r2, 02 , £2, 2:2) G K"* |a2 = -Aq} 
= {(r 2 , 02 ,£ 2 , 2 : 2 ) G K"* |a 2 = Ao} . 


Accordingly, we study the transition 


n 


- 1-62 

2 


A en,-|-b2 
^2 


^ A 


ex,-1-62 
2 


(^2)£2)^2) (i^2) ^2)-52)- 


For a qualitative description of X_^_i\^r 2 =o} and the objects defined above see hg. 11. 


(Ill) 


( 112 ) 



Figure 11 : Left: phase portrait of the corresponding layer equation of Center: phase portrait of the corresponding 

CDE of Right: Since the critical manifold is regular, by Fenichel theorywe know that the manifold r2 is perturbed 

to an invariant manifold r2,£2 which is at distance of order 0{£2) from r2. 

We know from section 2.1.2 that for sufficiently small £ 2 , there exists a change of coordinates that 
transforms X_^_i\[r 2 =o} into the vector field 


{ 02 = £2 

£'2 = 0 (113) 

Z!2 =-Z2, 

From the definition of the entry and exit sections (111), the time of integration is T = 2Ao. To obtain the 
component Z 2 of the transition need to integrate 

Z'2 = --Z2, (114) 

£2 

and then Z 2 = Z 2 {T). Therefore we have that after choosing a center manifold W'f, the transition H^^^ 
reads as 
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(115) 


n 7 f 
0 , £ 2,^2 exp 1-^—1 1 . 

Note that is an exponential type map. 

Analysis in the chart K_i 

In this chart the blow-up maps reads 

a = r|a 2 , b = —r^, 2 = r 2 Z 2 , s = r®£ 2 . 

Then we have that the blow-up vector field is given by 

-£2fr 

£ 2(1 + /aa) + £29a2 

£ 2/62 

— (z| — Z2 + 02) -|- £2/22 



nf^( 0 ,£ 2 ,Z 2 ) = 


(116) 


(117) 


where all the functions fe and ga 2 are flat along {r 2 = 0}. Observe that, as in the previous section, the 
subspace {r 2 = 0} is invariant. The restriction of X_i to this subspace reads as 


^—b\{r2=0} 


{ a'2 — £2 

£'2 =0 

4 =-{zl-Z2+a2). 


(118) 


The flow of X_i is a flat perturbation of the flow of X_i\^r 2 =o}- Therefore, let us continue our analysis 
restricted to the invariant space {r 2 = 0}. 


The manifold r 2 , which is defined by 


r2 = {(r'2,a2,£2,2:2) k2 = 0, £2 = 0, - ^2 -f 02 = O} 


is normally hyperbolic except at the two points p± = ± 


( 


2 

3V3’ 


1 

V3 


Let us define the sections 


A™’ = {(£2, 02, £2, 22 ) e I O 2 = -Ao} 

= I(r2,02, £2, Z2) € I 02 = 2I0} , 

where Aq > 0 is a sufficiently large constant. We are interested in the transition 

(£2, £2, 22) (£2, £2,52)- 


(119) 


( 120 ) 


( 121 ) 


For a qualitative description of X_i\^r 2 =o} and the objects defined above see fig. 12. 

Away from the fold points p±, the manifold r 2 is regular and thus, Fenichel’s theory applies. However, we 
need to take care of the transition near the fold point p+. The local transition of a slow-fast system near a 
fold point is investigated in e.g. [18]. However, in our current problem this transition is not essential. By 
this we mean that the passage through the fold point is seen as a flat perturbation of the trajectory along 
the stable branch of r 2 . In a qualitative sense, this is due to the fact that the transition H^^^ goes along a 
large NHIM, which fails to be normally hyperbolic only at one point. 
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Figure 12: Left: phase portrait of the corresponding layer equation of Center: phase portrait of the corresponding 

CDE of Right: The expected perturbed invariant manifold obtained from the flow of the corresponding CDE and 

layer equation. 


Proposition 3.3. We can choose appropriate coordinates in such that the transition II 2 : 

^en ,-62 , restricted to r 2 = 0 , is an exponential type map of the form 


U-’’%0,e2,Z2) 


^0, £ 2 , 4’2{s2) + Z 2 exp 



(Aq + £2'4’2{Z2, £ 2 )) 


( 122 ) 


where (j )2 are flat at £2 = 0 , ^/>2 is -admissible, and where ^0 is given by the slow divergence integral of 
^-bl{r2=0} ■ 

Proof. To prove that Aq is given by the slow divergence integral we proceed along the same reasoning as 
in proposition 3.2, so we do not repeat it here. In figure fig. 13 we see the three transitions that we must 
consider. 
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Figure 13: The three different transitions in which IIj is decomposed. The central transitions is locally an A 2 problem. The 
other two transitions at the sides are regular. 

The three transitions are defined as 



A®"-''2 ^ 17®" 


Y\fold 

^en ^ ^ex 

(123) 





where we define 17®" and 17®" as 
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( 124 ) 


— {(02, ^2: Z 2 ) G I 02 — —a2,en} 

= {(02, £2, Z2) G I ^2 = —^ 2 ,ex} , 

where 02 .en and Z 2 ,ex are sufficiently small positive constants. The total transition 11^^^ is given by = 
112*^®^ o o . Recall from appendix A that if we want to write the transition as an exponential 
type map, we require that is expressed as an exponential type map with no shift. The transition 112°*^ 

is studied in e.g. [13, 18]. In [13] is proved that there are local coordinates {Z 2 ,s) in 11®", and ( 0 - 2 , 5 ) in H®’', 
such that the transition is given by 


n{°*‘^(^2, £2) — (02, £2) 

— (^£2 + 0(62),£2J ■ 

Assume now that we have characterized an invariant manifold from H®" to fl®’' via the map 

Now we want to “extend” all the way up to the sections and via transitions along 

normally hyperbolic regions of r 2 . For this, it is more convenient to regard as a graph C 2 = (/)£ 2 (A 2 ) 

where (C 2 , A 2 ) are local coordinates around the fold point p+ and where (j)e 2 is a diffeomorphism for £2 > 0. 
In this way we can equivalently express the map 112°*'^ as 


n{°'"(C,£2) = (C2,£2) 

= (^’£2(0, £2) 


(126) 


where is a diffeomorphism for £2 > 0 and only a homeomorphism for £2 = 0. Next, following section 2.1.2 
we can find coordinates (Z 2 ,£ 2 ) in and coordinates (Z 2 ,£ 2 ) in in such a way that the 

transitions II^'^®’^ and II^'^®^ are given as 


nrn^2,£2) 


n”"®=(-Z2.ex,£2) 


Z2 exp-(Ag — 02 en) ) ) — (^ 2 , £2) 

' £2 

-Z2,exexp ^-^(Ag - 02 ) ) ) =(.^2,£2)- 


(127) 


Remark 3.4. Recall that along normally hyperbolic slow manifolds, it is possible to make a normal form 
transformation in such a way that this transformation respects certain constraint or structure of the vector 
field, [2, 3]. In this particular case, we respect the choice of the invariant manifold . 

Next, we can compute the composition II^o o by following appendix A and it thus 

follows that 


where 
at £2 = 0. 


n2-'’no,Z2,£2) 


^0, + Z 2 exp 


— (Ai + A 3 + £2'lp2^ 



(128) 


V>e2(0)exp(-f) 


and where 'ip 2 


fp 2 {Z,£ 2 ) is a C^-admissible function. 


Note that is Ant 


□ 
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3.6. Proof of theorem 3.1 

Let us first recall that, within the blow up space, we have three types of transitions according to the initial 
condition &i|a™, namely 

• If then we construct a transition passing through the charts iLgn Kg ^ Kg^. 

• If ^iIa^" ^ and ^iIa^" > 0 then we construct a transition passing through the charts Kg^ —>■ 

Kjgj) y Kg^. 

• If &i|a 5 " 4- and ^iIa^" < 0 then we construct a transition passing through the charts itTen —> 

K_-, ^ iLex. 

In fig. 14 we give a qualitative diagram of the local transitions obtained and their relationship. 

Let us only detail the transition through the inner layer A*”"®'' corresponding to &i|a'" € 0(^1^^)) th® other 
cases follow the same lines. 

The transition n'”''®'' : A*”"®'' —> Al’' is given as 

ninn®r ^ jjea ^ ^ex ^ ^ ^ jjinner (^29) 

where the matching maps are obtained from the blow-up map. For example, to obtain the matching map 
from the chart ATen to the chart Kg we relate the two directional blow-up maps 

a = —rf, b = rfbi, z = riZi, s = rlsi (130) 

and 

a = r|a 2 , b = r|52, z = r 2 Z 2 , £ = (131) 

Let us work out only with the z-component of the transitions as it is the only relevant one. Recall from 
section 3.2 that n™"®'' is an exponential type map with no shift. Next, the composition o 

Mgjj yields a diffeomorphism as II^^ is a diffeomorphism, and the matching maps are also diffeomorphisms 
on their domain of definition. Next, the last transition Ilg^ is an exponential type map with no shift, see 
section 3.4. Therefore, following appendix A we have that Ilg^ o ncentrai ^ jjmner jg exponential type map 
of the form 

n‘1"®’' = ^(Bi,£i) + Ziexp (A(i?i,ei) +ei^(i?i,ei,Zi))) , (132) 

where A > 0 and (p and are admissible functions. The differentiability of p and with respect 
to monomials is evident from the results of section 3.2. By blowing down we obtain that the transition 

Rinner . X;en ^ j^ex 

(in a small neighborhood of the cusp point and within the inner layer as domain) reads 
as 

n~ = <p{B,e) + Zexp(^-^{A{B,e)+e^{B,e,Z))^ . (133) 

To obtain the transition 11 : we now need to compose 11^“'' with exponential type maps on the 

left and on the right corresponding to 
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K--b 


Figure 14: All the transitions obtained in the charts. We have to compose all such transitions through the matching maps . 
A matching map relates the coordinates between the charts Ki and Kj. 
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n" : S- ^ E®" 
n+ : ^ E+. 


( 134 ) 


However, we must proceed with care. In order to express the transition H as an exponential type map, we 
need to choose appropriate coordinates on E“ and on S+ that respect the already chosen coordinates in 
S®" and in S®’'. Fortunately, this is possible with the extensions of Bonckaert [2, 3] to the normalization 
results of Takens [23]. 

For sake of clarity, let (Hen, ^en) be coordinates in E®" and (Hex, ^ex) be coordinates in E®’^. We have shown 
that these coordinates can be chosen in such a way that the “vertical” component of the transition map 

jjinner . y;en ^ 5.®^ reads aS 


n^en -^en, — -^ex 

= (()(Hen, 4” ^en exp ^ — ^^(Hen, s) 4“ (Hen, -^en)) ^ • 


(135) 


In this case the invariant manifold, say is given by Zen = 0. Using [2, 3] we can find suitable coordinates 
(H_, Z_) in E“ in such a way that 


n 2 _(H_,Z_,e) = Z_exp ^-^(Ho)^ = Zen- (136) 

In other words, there is a change of coordinates respecting the invariant manifold Alg under which the 
transition H” is an exponential type map with no shift and linear. Similar arguments hold for the choice of 

coordinates in E+. Finally, following appendix A, the composition H^^ o IIz^^ o leads to the result. 

A. Exponential type functions 

In this section, we discuss a particular type of function which will be found and used frequently throughout 
the main text. First, however, let us give two preliminary definitions. 

Definition A.l ( C^-admissible function). Let U € M". A function f : K" —)■ M zs said to be a C^-admissible 
function if f is -smooth away form the origin (for any £> 0), at the origin and if for all € N and 
Ui < I, there exists an N{ni) G N such that 

as 17,^0. (A.l) 

Now, we define a particular type of differentiability. For this we need to extend the common concept of 
monomial. In our context, a monomial, e.g. in two variables, u}{u,v) is any expression of the form u°‘v^ or 
of the form t6“(lnu)^, with a,/? € K. In general, if we let u G M™ and v G ffi", we allow a monomial u) to 
be any expression of the type ■u^’(lnu)'?, where and (Inu)"? = (Inui)'^i • • • (lnv„)*". We note 

that these monomials are admissible functions. 

Definition A. 2 (C^-function with respect to monomials). Let {U,V) G x M". IFe say that a function 
f{U,V) is a C^-function with respect to a monomial ui, if f is w.r.t. V in a neighborhood of 0 G M", 
and if there is a quadrant U = [0,mi) x ••• x [0,t6„) C M"* where the monomial uj is defined and such 
that the function f{uj,U,V) = f{U,V) is with respect to uj in 14. Similarly, the function f is said to be 
a -function with respect to the monomials wi,..., Wg if there is a quadrant 14 where the monomials are 
defined and such that the function f{uji,... , 0 Js,U,V) = f{U, V) is with respect to wi,..., Wg inlA. 
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Observe that a function / which is differentiable w.r.t monomials is an admissible function. As an example, 
consider f{U) = Ui\nUi(j){U) where (j){U) is smooth. This function is smooth away from U = 0 and (7° 
at the origin. However, it is not differentiable w.r.t. Ui at C/i = 0 but it is differentiable with respect to 
uj = UilnUi at w = 0. 

Let V G M™, Z G U, and as usual e denotes a small parameter. 

Definition A.3 (Exponential type function). A function D{V,Z,£) is called of exponential type if it has 
the following form 


D{V, Z,£) = B{V,e) + Zexp 




(A.2) 


where A and B, are admissible functions with A > 0, and B{V,0) = 0; and where is in z and 
w.r.t. monomials of {V,£) with ^(V,0,Z) = 0. We distinguish two particular cases 


1. The exponential type function D is without shift if B = 0. 

2. The exponential type function D is linear if ^{V, Z, e) = <f>(y,£). 


Remark A.l. Given a function D and if it is of exponential type, the representation of D is unique in the 
sense that all the functions in r.h.s of (A.2) are computable from D. In fact 


B = D{V,Q,e) 


A = lim 

Z-s-O 


$ = —e In 


f^ D{V,Z,e)-D{V,0,£) y^ 

DiV,Z,e)-D{V,0,e) \^ 


(A.3) 


We want to study the scenario where we have to compose D with some other functions and want to keep the 
exponential type structure. To be more precise, we consider D as an (V, e)—parameter family of functions 
(in Z) and compose it with a (E, e)—parameter family of diffeomorphisms on M. 

Proposition A.l (Composition on the left). Let : K —>■ K 6e a family of diffeomorphisms, and let 

D he an exponential type function. Then, the composition o D is also of exponential function of the 

form 


D = B{V, e) + Z exp 


A{V,£) + ^iV,Z,e) ^ 


(A.4) 


where B and are admissible functions. 

Proof. Let us simplify the notation by writing = ^'(ye). Since 4* is a diffeomorphism we can write 
'I'(a + b) = 'l'(a) + (7(1 + V'(a, b))b, near 6 = 0, with ip a function such that ip{a, 0) = 0 and with (7 > 0. 
Then we have 


d;* o D{z) = 4/ 


^B + Z exp 


7l + $ 


= 'S(B) + (7(1 + ip(V, Z,£))Z exp 




(A.5) 
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Since C > 0 we can take the logarithm of (7(1 + Z, e)) and then we have 


d/ o D{z) 


d>(S) + exp(ln((7(l + '4^))Z exp 




'SiB) + Zexp 


^ + $ + eln((7(l + V’) ^ 


(A.6) 


The result is obtained by setting B = 'i’{B) and $ = <& + eln((7(l + tp). □ 

Proposition A.2 (Composition on the right). Let ^'(v,£) : M —)■ M be a family of diffeomorphisms with no 
shift, that is 'l'(y_e)(0) = 0 for all (V)e), and let D be an exponential type function. Then, the composition 
D o is also of exponential function of the form 


D = B{V, e) + Z exp 


A{V,e) + ^iV,Z,e) ^ 


(A.7) 


where B and are admissible functions. 

Proof. Let us simplify the notation by writing Since 'I'(O) = 0 we can write 'I'(z) = C(l + 0(z))z 

with C > 0. Then we have 


Do'^{z)= D{Cil + 0{z))z) = B{V, e) + (7(1 + 0{z))z exp - 


AiV,s) + ^V,e,-^) 


= B{V,e) + 2 :exp - 


A{V,e) + ^{V,e, '!') + eln((7(l + 0{z))) 


(A.8) 


The result then is obtained by setting $ = <i)(P, e, 'L) + e ln((7(l + 0{z))). □ 

Remark A.2. If we want the composition 11 o to be of exponential type, the family cannot be 

arbitrary. In order to preserve the “exponential structure”, 'I'(y,£) should satisfy the hypothesis of proposi¬ 
tion A.2. In corollary A.2 we show a particular case in which the diffeomorphism 'k can have a shift and 
yet preserve the structure of the exponential type function. 

Let us proceed by presenting a couple of useful corollaries. 

Corollary A.l. Let Di and D 2 be two exponential type functions of the form 


Di (V, Z,£) = Z exp ( — 


AiiV,s) + ^iiV,Z,£) 


D2{V,Z,£)=B2{V,e) + Zexp - 


A2{V,e) + <I>2{V,Z,£) 


(A.9) 


that is, Di is an exponential type function with no shift. Then D 2 o Di is an exponential type function. 
Corollary A.2. Let Di and D 2 be two exponential type functions with D 2 linear, this is 

Ai{v,£) + <i>i{v,z,£y 


Di{V, Z, £)=Bi{V,£) + Z exp 1^- 
D 2 (V, Z, £)=B 2 iV,£) + Z exp (- 


A2iV,£) 


(A.IO) 


Then the composition D 2 o Di is of exponential type. 
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It is useful to consider the following: let X{y,Z^e) be a given vector field on and let Sq and Si 

be codimension 1 subsets of which are transversal to the flow of X. For the moment it is sufficient 

to think of a section given by {Vj = vq} or by {e = Eq} with vq and Eq fixed constants. Induced from 
definition A.3 we then have the following. 

Definition A.4 (Exponential type transition). A transition 11 : Eq —>■ Si is called of exponential type if 
and only if its Z-component is an exponential type function. This is, an exponential type transition is of the 
form 


IliV,Z,e) 


{G,D,H) 

(G(V,e), B{V,e) + Zexp 


AiV,e) + ^V,Z,e) \^ 


H{V,s) 


(A.ll) 


where G : —>• and H : —>• M are with G(V, 0) = E and H{V, 0) = 0; where A, B and $ 

are -admissible functions. The names exponential type transition with no shift and linear are inherited as 
well from the type of D. 

Suppose now that A is a given vector field on as above, and let S^ with i = 0,1, 2, 3,4, 5, be disjoint 

sections which are all transversal to the flow of X. Assume that X induces exponential type transitions 
Hi : Si_i —>■ Si with i = 1, 2,3,4, 5 of the following form 


1 . Hi is with no shift and linear 

2 . 112 is with no shift 

3. Ha is a general diffeomorphism 

4. 114 is with no shift 

5. Ha is with no shift and linear. 


We need to show that the composition of all these five maps is an exponential type transition. 


Proposition A.3. Let Hi : Si_i —>■ S^ as described above. Then the composition = 1150114011301120111 
is an exponential type map of the form 


n = 


^G(14,e), B{V,e) + Zexp 


A{V,s) + ^{V,Z,e) ^ 


H{V,e) , 


(A.12) 


where A = Ai -\- A 2 + A 4 + A 5 . 

Proof. Let us write each of the transitions as follows. 

1. ni(E,Z,E) = (Gi,Di,Hi) = (Gi,Zexp(-^^i^),iLi) 

2. U2{V,Z,e) = (G2,D2,H2) = (G2,Zexp 

3. n3(E,Z,E) = (G3,D3,H3) 

4. n 4 {V,Z,s) = (G 4 ,D 4 ,g 4 ) = (G 4 ,Zexp (^_ AAV,e)+^AV,Z,e) ^ 

5. n^{V,Z,e) = {G^,D^,H^) = (G^,Ze^p(-^^),H^) 
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For brevity let 112 o Hi = (G 2 , D 2 , i? 2 )- Then we have 


(G 2 , £* 2 ,-^ 2 ) — 

(^G2(Gi,i/i),Diexp 


( ^2(Gi, i7i) + $2 (Gi,D i,iJi)\ 


,i?2(Gi,i7i) . 


(A.13) 


Now, we take care only of the Z-component of the composition 112 o IIi. From the hypothesis on Gi and 
Hi we can write Gi = V + 0{e) and Hi = Q;e(l + 0(e)) with a > 0, then 


D 2 


Z exp 


^ Zli(V,s) + A 2 (V,e) + ^ 2 iV,Z,s)'^ 


(A.14) 


where we have gathered in <I >2 the function <I>i and the terms resulting from taking Gi = V + 0{e) and 
Hi = ae{l + 0(e)). In a similar way, letting IIs o 114 = (G 5 , we get 


G 5 


Z exp 


^ ^ 4 ( 5 ) +-45(g) + $ 5 ( 1 ^, Z,g) ^ 


(A.15) 


Next, and following similar arguments as above, we know from proposition A.l that the composition 11321 = 
Ha o 112 o III is of exponential type with shift. Finally since the transition 1154 = Ha o 114 is of exponential 
type with no shift, and using proposition A.l, we have that 1154 o 11321 is an exponential type transition as 
claimed in the proposition. 

Remark A.3. In the case where Ha is an exponential type map, we get a similar result with A = Ai + 
A2 + As + A4 + As ■ 


□ 


B. First order differential equations (by R. Roussarie) 

The contents of this section shall appear in greater detail in [ 6 ]. We reproduce some results here for 
completeness purposes and to use them in appendix C.l. 

Let X{x) be a smooth vector field defined on IF C K”, for arbitrary n S N (here we include the possible 
parameters). Let G{x, y) : IF x K —)■ M be a smooth function. We shall study the solutions of the first order 
differential equation 


X ■ K{x) = G{x,K{x)), (B.l) 

where K(x) is the unknown function. We assume the following 

1. There exists an open section S C IF which is transverse to X. 

2. Let (f>{t,x) denote the flow of X. We can choose an open domain IFs with the property that for any 
X € IFs, there exists a unique smooth time t(x) (possibly unbounded) such that (j){t{x),x) G S. 

3. The vector held Z{x,y) = X{x) + G{x,y)dy has a complete flow. 

The flow of Z takes the form {(j){t, x),'tf{t, x, y)), where (f) is the flow of X. It follows that K{x) is a solution 
of (B.l) if and only if the graph {y = A'(x)} is a surface tangent to the vector field Z. Then we have the 
implicit formula 


tp{t{x),x, K{x)) = 0. 


(B.2) 
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In our applications, the function G is affine in y, that is G{x,y) = L(x)y where L and 11 are smooth. 

If we write L{t,x) = L{(j){t,x)) and fl(t, a;) = Il{(j){t,x)) (where (j) is the flow of X), we have for ip the 
following linear differential equation 


dpj 

dt 


{t, X, y) = L{t, x)'tp{t, X, y) + n(t, x). 


(B.3) 


Then we can integrate (B.3) with the initial condition 'ip(Q,x^y) = y to obtain 


= exp ( ^ L{T,x)dT^ ^ J 


exp ( — / L(f 7 , x)da 

Jo 


dr > . 


(B.4) 


Since exp ^ J* L{T,x)dT^ > 0 we can solve the implicit equation (B.2) obtaining 

ft{x) 

lO I \ Jo 


ftG) [ / 

K{x) = — J n((/)(r, a;)) exp i—J L{(j){a, x))da 


dr, 


(B.5) 


where we recall that (p is the flow of X and t{x) is the time to go from x to the section E along this flow. 

Let us now assume that the vector held X is partially hyperbolically attracting in the following sense: we 
assume coordinates a; = (a, 6) and that the vector held X has a decomposition X(a;) = U{x)+V{x) 

where U is the component along and V is the component along M"?. Moreover, we assume that V = Q 
on RP X {0} (that is X is tangent to Rp x {0}). We also assume that at each point x = {a, b) it is satished 
that DbV{a, 0) has all its eigenvalues with strictly negative real part. We further suppose that X is given on 
W = D X A where D is a domain diffeomorphic to a ball in and A is a ball in We choose A = Ap^ 
for some po > 0 where Ap = {5 G | 11&| | < p}. It then follows that under a linear change of coordinates 
(a, 6) i-G {a,A{a)b), the vector held X enters along D x dAp for 0 < p < po if we choose po small enough. 
We now have the following 

Proposition B.l. Assume that DbV{a,0) has all its eigenvalues with a strictly negative real part and that 
Po is small enough as explained above. Let B be any domain diffeomorphic to a closed ball inside the interior 
of D and assume that the function n(a;) is flat along D x {0}. Then the equation 


X ■ K{x) = L{x)K{x) + n(a;) 


(B. 6 ) 


has a smooth solution K(x) in B x A which is flat along B x {0}. 

Proof. Let f{a) : Rp —>■ [0,1] be a smooth function which is equal to I on B and equal to 0 on a neighborhood 
of dD. Dehne the vector held 


T = V + fU. (B.7) 

This vector held T coincides with X on B x A. Moreover, T is tangent along dD x A and enters the domain 
D X A along D x dA. Let 4>{t,x) = {(j)a{t,x),(l>b{t,x)) G x denote the how of T. It follows that 
(j){t, x) G D X A for all a; G D X A and all t > 0. From the assumption on V we have that there exists a 
positive constant E > 0 such that 


\\M't,J:)\\ < \\b\\exp{-Et), 


(B. 8 ) 


for any x = {a,b) G D x A and t G [0,+oo). We now want to use this how (p in (B.5) noting that if the 
integral converges, then K{x) is a solution to the equation T ■ K = LK + 11 on I? x A and then to the 
equation X ■ K = Lk + 11 on B x A. In this setting (B.5) is written as 


K{x) 


n(())(T,a;)) 


exp - 


L(^(a, xy)da 


dr. 


(B.9) 
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Now, we need to prove that (B.9) defines a smooth function on Z? x A which is flat along D x {0}. In other 
words, we shall prove that K and all its partial derivatives are equal to 0 on Z? x {0}. As L is bounded, 
there exists a constant Mq > 0 such that 


exp 


L{(l>(a, x))da 


< exp(Mor). 


Next, let iV € N. Since n is flat in v, there exists a constant Pn > 0 such that 




(B.IO) 


(B.ll) 


and then from (B.8) it follows that 


|n((/)(T,a;))| < PA,||6||^exp(-AAr). 


(B.12) 


Using these estimates we have that 


r+ca 

\K{x)\ < PAr||6|r / exp((Mo - NE)T)dT. (B.13) 

Jo 

The integral in (B.13) converges if N is large enough, strictly speaking if A > This proves that by 
choosing N sufficiently large, the right hand side of (B.9) defines a function which is continuous and equal 
to 0 on I? X {0}. 

Let us now consider any partial derivation daK of K. Let us write 


H{t, x) = x)) exp 


L{(l>(a, x))da 


(B.14) 


the integrand in (B.9). Using chain rule on the derivative of (B.9), we have to prove that the integral 


+ 00 

daH{T,x)dT (B.15) 

is convergent and that there is an estimate similar to (B.13) for N large enough. We do not want to give all 
the details here and refer the reader to [6]. The idea is that daH{T, x) is a finite sum of terms such that each 
of these terms is a product of factors which are partial derivatives in x and are of of one of the following 
forms 


1. dai{<p{T,x)). Since 11 is smooth and fiat along D x {0}, this is also the case for dai{4>{T,x)). Therefore, 
for N sufficiently large, we can write an estimate of the form 

\daAHT,x))\ < PN^J\b\\^ expi-NEr), (B.16) 

for constants Pn^^ > 0. 

2. da24‘{T,x) (resp. 9 q,2^((t, x), note that 0 < ct < t). By the usual variational method along trajectories, 
there exists constants E^^ > 0 such that \da 24 ’{x,x)\ < exp{Ea 2 T) (resp. |i9a2^(CT,x)| < exp(Uc, 2 (T) ). 

3. 9 q3L(^(t, x)). As L is smooth in I? x A, all these factors are bounded by a constant Mq-j. 

4. exp (— fj L(^(i7, x))do'). This factor is bounded by exp(MoT). 

Next, by remarking that a factor of the first type appears in each term of the expansion of daH, and taking 
N large enough, it is possible to conclude that the integral (B.15) converges an is equal to 0 for x € I? x {0}. 
Therefore, the partial derivative daK{x) exists, is continuous and is equal to 0 on P x {0}. 

□ 
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C. Normal form and transition of a semi-hyperbolic vector field 

In this section, we present a rather general framework for the computation of a normal form and the 
corresponding transition of a vector fields with a semi-hyperbolic singularity. The contents of this section 
are not only relevant for the object studied in this document, but for more general systems as well, c.f. [13]. 
To make our computations simpler, we prove a lemma that allows us to “partition” a smooth function. As 
a simple example of this partition, let f{u,v) be a smooth function on We show that / can be written 
as f{u,v) = fi{uv,u) + f 2 {uv,v), where /i and /2 are smooth. This type of result becomes useful when 
computing the transition map that we present in appendix C.3. 

C.l. Normal form 

Here we provide a normal form of a semi-hyperbolic vector field which frequently appears in the analysis 
of slow-fast systems. The goal of obtaining such a normal form is that the computation of the corresponding 
transition becomes simpler. 

Proposition C.l. Let a, j5 = {j3i,. .., Pm) and 7 be non-zero constants, and consider the vector field X 
given by 


X : 


z 


/ 


= awu{l + /) + wg 
= PjWVj{l-\- f) 

= 7u; 2(1 -h /) 

= —A -|- h, 


(C.l) 


where j = 1 , 2 ,..., to; where the functions f = f{u, v,w,z), g = g{u, v, w, z) and h = h(u, v, w, z) are smooth 
functions which are fiat at the origin and where A = K(u,v,w,z) is a smooth function such that 

A(0) = 0 and § 7 ( 0 ) > 0. Then there exist a coordinates {U,Vi,... ,Vm,W, Z) under which X can be 
written as 


U' 

= aWU 

V' 

3 

= P3WV, 

W' 

= 71^2 

Z' 

= -GZ, 


(C.2) 


where G = G{U, V, W) is a G^ function such that G(0) > 0. 

Proof of proposition C.l. From the definition of the vector field X we note that the origin is a semi- 
hyperbolic singular point. The hyperbolic eigenspace is 1-dimensional while the center eigenspace is (to- 1-2)- 
dimensional. We now proceed in 4 steps as follows. 

1. Define a new vector field Y by Y = j^X, which reads as 


{ u' = awu wg 

vi = PjWVj 

w' = 

z' = —A -I- h, 


(C.3) 


where the functions g and h are flat at the origin of Note that in a small neighborhood of 

{u,v,w,z) = (0,0,0,0) the vector fields X and Y are smoothly equivalent. 
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2. By looking at DY{0), there exists an (to + 2)-dimensional center manifold [11]- Let Mq be the set 
of critical points of Y, that is 


Mo = {(u, V, w, z) I A{u, V, 0, z) = 0} . (C.4) 

By definition, the manifold Mq is invariant and normally hyperbolic. Now, assume |w| <C 1. This 
condition appears naturally in our applications. By Fenichel’s theory [10] the manifold Mq persists as an 
invariant normally hyperbolic manifold M^,, for sufficiently small w 0. We identify with yV ’2 ■ In. 
other words, there exists a function to = m{u, v, w) such that the center manifold W 2 is given as a 
graph 


yy 2 = Graph(M, u, w, to). (C.5) 

Define = z — m, then C^' = z' — m!. But we know, due to invariance of under the flow of Y, that 
C'|(=o = 0. This is, there exists a function H = H{u,v,w,Q such that C' = —H(. With H{0) = 0 

and ^(0) > 0. 

dC 

In conclusion of this step, there exists a transformation ijj : {u,v,w,z) 1 —>■ {u,v,w,Q that transforms 
the vector field Y into 


{ u' = awu + wg 

= PjWVj 

w' = 

C =-HC, 


(C.6) 


OH OA. 

where H = H(u,v,w,Q is a function such that iL(0) = 0 and where (0) = v^(0) > 0. 

az az 

3. Observe that thanks to the previous step, the center manifold has the simple expression = 
{C = 0}. We now want to separate the variables on the center manifold (these are {u,v,w)) from those 
on the hyperbolic subspace (z). Additionally, we want to keep the simple format that Y has in the center 
direction. This amounts to find a change of coordinates along C only. For this we use an extension of 
Takens’s theorem on semi-hyperbolic vector fields [23] due to Bonckaert [2, 3]. With this, it is possible 
to show there exists a transformation, fixing the center coordinates, that conjugates Y to the vector 
field 


{ v! = awu + wg 
v'. = /3jWVj 

w' = 

Z' = -HZ, 


(C.7) 


where now the flat perturbation g is independent of Z and H = H(u,v,w) is a function with 
H{0,0,0) > 0. 

4. In this last step we eliminate the flat perturbation from Y, which appears only along u. Due to the 
previous step, the dynamics on the center manifold are independent of Z. The restriction of Y to W 2 
reads as 
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u 


' = awu + wg 

Note that for w 0, the vector field ^Yly^c is hyperbolic. Let y = ^Y\y^c, that is 


(C.8) 


{ u' = au + g 

V'j = PjVj 

w' = JW. 


(C.9) 


Now we have the a result that shows that there exists a change of coordinates, respecting the variables 
{v,w) that kills the term g. Keeping the coordinate w fixed is important because we want to prove an 
equivalence relation with wy and not with y. The following proposition shall appear in a general context 
in [6]. 

Proposition C.2 ([6]). There exists a diffeomorphism {u,v,w) i—)■ {u + H{u,v,w),v,w) with H flat at 
(u, v,w) = 0 which brings y to 


u' 

= an 



II 

(C.IO) 

w' 

= 7'u;. 



Proof. We shall use the path method to show that y is conjugate to y. Let s be a parameter and let us 
define the s-parameter family of vector fields 




y + sg 


d_ 

du 


(C.ll) 


We call jV® the path between y and 3^ + We now look for an s-parameter family of diffeomorphisms 
TL^ with TL'^ = Id such that for each s we have the conjugacy 

Hiy = y^. (C.12) 


In such a case, the vector fields y and y + g^ are conjugated by . By derivation of the family "H® 
along s, we obtain an s-parameter family of vector field C® satisfying 


C®(H®) 


aiL® 

~d7' 


(C.13) 


This implies that by derivation of (C.ll) with respect to s we obtain 


Therefore, if are able to find a solution C® of (C.14), the conjugacy H® is obtained by integration of 
(C.13). In our particular case, we are looking for a solution along the u-direction, that is of the form 
d 

C® = Ps 7 ^ ■ It follows that 
ou 
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[r,ci 


(C.15) 


9 d d 



' du 



JL 

du' 


Therefore we have reduced our conjugacy problem to solving the differential equation 


y^{Ps) - Ps = 9, (C.16) 

where we recall that g = g{u,v,w) is flat at {u,v,w) = (0,0,0). We now want to use proposition B.l 
to show that (C.16) has a solution Pg = Ps{u,v,w) which is flat at {u,v,w) = (0,0,0). For this, let 
Gs = a + sff • Now, we only need a small adaptation: in the setting and notation of proposition B.l 
we may assume (under the suitable arrangement of coordinates) that (or X in proposition B.l) is 
tangent to x {0} and {0} x Let A4^(a) and A4^(b) denote the space of germs of s-families of 

smooth functions that are flat at {a = 0} and at {& = 0} respectively. Using a blowing-up at 0 € M" it 
can be shown that A^“(a, b) = A4^(a) + A4^(b) (see the arguments in lemma C.l). From this formula, 
it follows that it is sufficient to solve (C.16) in the spaces M^(a) and Mf{b) respectively. Naturally, 
these two cases are equivalent up to the change of by —3^® and Gg by —Gg in (B. 6 ). In either case, the 
vector field (or —3^®) of (C.16) satisfies the hypothesis of proposition B.l. Then for g in Mf{a) (resp. 
in A4^(b)) and applying proposition B.l, we can solve (C.16) with Pg in A4^(a) (resp. in M^{b)). 

□ 

Thus, from proposition C.2, we have that 3^ ~ 3^ respecting w, which implies wy ~ wy. Therefore, we 
conclude that (C.7) can be written as stated in the proposition. 

□ 


C.2. Partition of a smooth function 

In this section we investigate the problem of partitioning a smooth function. The result presented below is 
important since it is used to simplify the computation of transition maps. To be more specific, let us give a 
brief example. Consider the three dimensional differential equation 

x' = X 

y' = -V (C.17) 

z' = 9{x,y)z, 

where g is a smooth function. We want to take advantage from the fact that xy is a first integral. We show 
below that the function g can be partitioned as g(x,y) = gi{xy,x) + g2ixy,y). This makes the integration 
of z' simpler. 

Lemma C.l. Let m € M and v € M™. Let f = f{u,v) be a smooth function such that /(0,0) = 0. Then 
there exist smooth functions fo = fo{uv,u) and fi{uv,v) such that the function f can be written as 


f — fo + fi, 


(C.18) 


where /o(0, 0) = 0 and /i(0, 0) = 0. 

Proof of lemma C.l. We proceed in two steps. The first consists in proving the formal version of the 
statement. The second step is to extend the formal result to the smooth case. 
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Formal step 

Let / denote the formal expansion of the smooth function /. Let p € N and q G N™. We use the following 
notation: 

• By 9 > 0 we mean Qi > 0 for all i G [1 ,to]- 

• For a vector v G M™ we write v'^ = . 

• The Li norm of q is denote by |< 7 |, and thus for g > 0 we have jg] = J^JLi Qj- 

• We denote by qi the vector 

Qi — {Qlt ■■■ J Qi—lt Qi+lt ■ ■ ■ ! Qm) (C.19) 

and therefore we have that v^' reads as 


— 


^qi+1 


Besides, we have that the Li norm of qi is given by \qi \ = jg] — qi = 




(C.20) 


The formal series expansion of / reads as 


/ - ^ OpgU^V'^, (C.21) 

p>0,g>0 

where ooo = 0. With the notation introduced above, we can partition / as follows 

m 

/ = H H apg{uVi)P{v,vy'v‘l'~^~^'^'^, (C.22) 

p>\q\ ®=i 9 i>p+l 9 d 

where 


{uvf = ■ ■ ■ {uv^Y- 

{v,vY' = 


(C.23) 


and where a'pg G K are suitable chosen coefficients. Let r G N™, s G N. Define the following formal 
polynomials 


h{uv,u) = ^ ars{uvYu’^ = ^ a'pg{uvYu^ 

r,s>0 P^\q\ 


where € M, and 


gi{uvi,v) = ^ Pirsiuv^Yv'' 

r,s,i>0 

= a!pg{vivY'{'U‘ViYvf~^~'''^''' 

9i>p+|9i| 


(C.24) 


(C.25) 
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where j5irs € IR. The coefficients ars and j5irs are conveniently chosen to make the definitions hold. Let 
uv = {uvi,..., uvm)- Define g = g{uv, v) by g{uv, v) = v), then we can write / as 


f{u,v) = h{uv,u)+g{uv,v). (C.26) 

This shows that the proposition holds for formal series. 

Smooth step 

By Borel’s lemma [5], there exist smooth functions h = h(uv,u) and g = g{uv,v) (whose formal series 
expansions are h and g respectively) such that 


f = h + g + R, (C.27) 

where R (reminder) is a flat function. We now show the following. 

Proposition C.3. Let u € R, v € M'", and R{u,v) be a smooth flat function at (0,0) G M x M™. There 
exist flat functions tq = ro{uv, u) and ri = rfluv, v) such that 


R = ro +ri. 


(C.28) 


Remark C.l. Proposition C.3 together with the formal step f = h + g imply our result. 

Proof of proposition C.3. For this proof we shall use the blow-up technique. Let $ : S'™ X M+ ^ 
be a blow-up map. The map <I> maps S™ x {0} to the origin in Let i? be a function defined by 

R = R o Since R is fiat at the origin, the function R is flat along the sphere S™. We assume that the 
function R = R{u,v) is defined on a small neighborhood TZ of the origin in M x M™; this neighborhood is 
defined as 


n = {\u\ < A, \v^\ < B,} , (C.29) 

for some A,Bi positive scalars. Let 0 < d < 1. The sphere S™ can be partitioned into m -I- f regions as 
follows: 


7/ = S™\{|m| <d} 
V. = S™\{|t;,| <d}. 


(C.30) 


where {u, v) = (u,vi,..., Vm) € S™. We can then take a partition of unity to split R as 


R{u,v) = Ro{u,v)+'^Ri{u,v), (C.31) 

i=l 

where Supp(i?o) C U and Supp(i?i) C Vi for i G [1,to]. We define as Rq and Ri the corresponding functions 
on flat at the origin given by the blow-up map <&, that is Rj = Rj o <I), for j = 0,1,..., m. Note that 

i? —> i? is an isomorphism between the space of functions on {u, v) G flat at the origin, and the space 

of functions on {{u,v),p) G 5"™ x IR+ fiat at 5"" x {0}. Therefore, the splitting (C.31) induces the splitting 


m 

R{u,v) = Ro{u,v)+ '^Ri{u,v) (C.32) 
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(C.33) 


of functions on We will now prove that there exist flat functions tq and such that 

Roiu,v) = roiuv,u) 

Ri{u,v) = ri{uvi,v). 

Let us detail only the case of Rq. The other functions are obtained in a similar way. 

The function Rq has support in U. We can parametrize U by the directional blow-up map which reads 
as 


{u,Vi,..., Vm) (m, MWi, . . . , UVm) = (m, Ui, . . . , Vra)- (C.34) 

Now, suppose that there exists a flat function Pq defined by 


Ro{u,v) = Po{u,u^v). (C.35) 

This implies that there is a function tq = Pq o <f>“i such that 

Ro{u,v) = ro{u,uv), (C.36) 

which is precisely what we want to prove. So, now we only need to show that indeed a function Pq as above 
exists. For this let us define coordinates {U,Vi,..., Vm) given by 


U = u,Vi= u^vi, ...,Vm = v?Vm, (C.37) 

and let Po(u, V) be a function defined as 

Pq{u,V)=Rq[^,^. (C.38) 

Note that Pq is flat at {u,V) = 0. This is seen as follows. Since Rq is flat along {u = 0}, it follows that 
Po(0,0) = Po|„^o = 0 and 


dpQ 

du 


( 0 ) 


dpQ 

dV 


( 0 ) 


u=0 


dRQ 
du 
1 ORq 
dvi 


= 0 


Lt—0 — 0; 


(C.39) 


and so on for the higher order derivatives. 

Finally, for convenience of notation we define rQ{uv,u) = rQ(u,uv), thus we can write Rq{u,v) = rQ{uv,u) 
Following similar arguments as above we find the functions = ri(uvi,v) such that Ri(u,v) = ri(uvi,v) for 
i G Then we define ri{uv,v) = YPiPi ''’i{uvi,v). It follows that 


R{u, v) = tq^uv, u) + ri(uv, v). 


(C.40) 


□ 


With this last proposition we can now write the function / as 

/ = h{uv, u) + g{uv, v) + R{u, v) 

= h{uv, u) + g{uv, v) + ro(uw, u) -I- ri{uv, v). 

Finally, to show the lemma we define the smooth functions fi, fi of the statement by 

/i = -P ro 
]’i=gPri. 


(C.41) 


(C.42) 

□ 
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C.3. Transition 

In this section we investigate the transitions for the vector field computed in appendix C.l. Relabeling 
the coordinates we recall that reads as 


u' 

= awu 


= jdjwv 

w' 

= 'yw'^ 

Z' 

= -gZ, 


(C.43) 


where j = 1,2,... ,m, and where g = g{u, v, w) is a function such that (/(O) = A > 0. We assume that 

w € M“'". For our applications, we are interested in only two particular situations. 

1. The saddle 1 case where a = — 1, > 0 for all j € [1, to], and 7 > 0. 

2. The saddle 2 case where a = 1, j3j <Q for all j G [1, to], and 7 < 0. 

Saddle 1 

In this case we investigate the transitions of a vector field of the form 


U 



U' 


= —wu 

= PjWVj 

= 'yw'^ 

= -9Z, 


(C.44) 


where the coefficients j3j, 7 are positive. Observe that the flow in the direction of u and Z is a contraction 
while it expands in all the other directions. Roughly speaking, this implies that a transition can go out at 
any expanding direction vj of w. 


We investigate two types of transitions that are used in our applications. For this, let us define the following 
sections 


Sen = {(m, V,ia,Z)lu = uj 
= {(^t,V,W, Z)\w = Wout} 

U, W, Z) 1 Vj = Vj^out} ■ 

In this section we compute the transitions 

n"' :Se„ ^ S“ 

(v,w,Z) {u,Vi,Z), 


for all i G [i, to], and 




:Si —)■ S 


±Vj 

ex 


(v,w,Z) {u,Vi,W,Z), 


(C.45) 


(C.46) 


(C.47) 


for all i G [I, to] with i ^ j. 

Proposition C.4. Consider the vector field Y given by (C.44) and let Sen, S^“^ and II™, 11=*=*'^ he as 
above. Then 
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• The transition 11“ is given by 


w \ 


1/7 


u = u ^ 

^ '^out 

Z = Z exp 


; ■ 
(> 


Vj = V. 


'^out \ 


t Wout \ 
\ W / 


A 


-(1 + aw ln(w) + wG 

jw ' 


(0.48) 


where a = a{uvy^%uw^'^'^) and G = G(uvl^^\uw^^~*, yii) are G^ functions with Hi = v^^'w . 

The transition is given by 


u = \ — 

.Vj 


1//3 


Vj 

Vi=V^\ — 


ft/ft 


w = w \ — 

.Vjj 


vA 


7/ft- 


Z = Z exp 

with i ^ j and where 


A 

jw 


^1 + a'w In(wj) + wG' 


(0.49) 


a' = 01 {uvA' ,uw^^^) 

G' = G'{uvA^fuw^^'*, hi) 


(0.50) 


are functions with hw = and hi = ft • 

Proof of proposition C.4. We detail first the computations for the transition 11“. The transition 
is computed in a similar "way so we only highlight the key parts of the computation. 

The transition 11“ 

In this case, the time of integration is T = In , where Wout = 'M^(0lsft ^ “ ii'(^)|s,„- This time 

of integration is obtained form the equation w' = jw. We also make the assumption that Vi € 0(w^’^'^). 
This assumption appears our applications, but roughly speaking it ensures that Vi is well defined when 
w —)■ 0. From the form of Y we evidently have 


u(T) = u = u 
Uj(T) =Vi = v, 


( "" 1 

1/7 

\Wout) 


(Wont 

\ ft/7 

H ( 

\ w ■ 

) 


(C.51) 


It only remains to compute the transition for the Z coordinate. Let us rewrite Y as follows 

u' = —u 

ft = Aft 

w = yw 

^ KyG{u,v,w) ^ 
w ’ 


(C.52) 


where G is a function vanishing at the origin. Observe that we have the first integrals u^'Vi and u'^w. 
We shall take advantage of such a fact. We define new coordinates {U, V, W) given by 


U = u, Vf' = Vi, = w. 
In these new coordinates we have the system 


(C.53) 
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(C.54) 


U' = -U 
VI = V 
w = w 

VFT- 

In the new coordinates, the time of integration is given as T = In (^)- To have an idea of the expression 
of Z, let us first study a simplified scenario. 

The case G = 0 

Let us suppose G = 0. Therefore we have Z' = which has the solution 


Z(t) = Z( 0 ) exp A J W{s) '*ds^ , 


(C.55) 


where W{s) = IT(0) exp(s). Substituting the time of integration T we have 


Z(r) = Z = Zexp - — 


A 


A 


= Z exp — 


7IT')' 


1 - 


e 

) 

w: 


(C.56) 


Observe that Z —>■ 0 as W 0. Let us now study the general case. We expect that the general case G 7 ^ 0 
is a perturbation of (C.56). 

The case G 7 ^ 0 

We now consider that G 7 ^ 0, we have 


Z(T) = Z = Zexp (lo + Ii), 


(C.57) 


where 


/n = -A 


]V(s) 


ds 


h = 


G{U{s),V{sf^,W{sV) 

W(s)T' 


(C.58) 


ds. 


The integral /q has already been computed above. Let us write F{U,V,W) = W(^) 

do this because G(C/, 0, 0) = 0 and e 0{W'^). Now we estimate the integral /i. Using lemma C.l, we 
can write 


h = [ [Fi{s) + F 2 {s)]ds, 

JO 


(C.59) 


where 


Fi = FiiUVi, UVm, UW, U) 

F 2 = F^iUVi, UVm, UW, Ui, ..., Vm, fU). 


(C.60) 
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Observe that UW and all the UVj's’ are first integrals. Let Ji = J Fi and J 2 = J F 2 - Then we have 

Ji = / Fi{UV,UW,U{s))ds 
Jo 

= / Fi{UV,UW,Ue-^)ds. 

Jo 

Let us make the change of variables y = e~‘, we obtain 

Ji = - Fi{UV,UW,Uy)^. 

Ji y 

We expand the function Fi in power of y that is 

Fi{UV, UW, Uy) = Fi{UV, UW, 0) + 0{y). 

Then we have 

W 

/• vvo dy ~ 

Ji — — ai - \- Fi, 

Ji y 


{CM) 


(C.62) 


(C.63) 


(C.64) 


where ai = ai{UV, UW) and Fi = Fi{UV, UW, Uy{T)) is some (unknown) function. Finally we get 


ITo 


W 


J, = a,ln ^ +F, UV,UW,U— . 


W 


ITn 


(C.65) 


The function Fi is but unknown, and Wq is a fixed positive constant, then we can simplify the notation 
of Fi as Fi =Fi{UV,UW). 


Next we have 


J 2 = f F 2 {UV,UW,V{s),W{s))ds 
Jo 


= / F2{UV,UW,Vie^^\...,Vme^^^We^^)ds. 

Jo 

Let us make the change of variables y = . Then we obtain 

Wq 

J2= f'^ F2{UV, UW, Fi/i,..., F™/™, Wy'')^. 

Ji y 


As above, we expand in powers of y, that is 


(C.66) 


(C.67) 


F 2 = a 2 + 0{y), 


(C.68) 


and then we have 

J2=a2\n(^^^+F2, {cm) 

where q ;2 = a 2 {UV,UW), F 2 = F 2 {UV,UW, yi) is a function with yn = ViW~^ for all i G [l,rn]. Recall 
that since Vi G 0{w^'/^) we also have that V G 0{W), that is fXi is well dehned. 
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Now we can write the integral Ii as 


/i — Ji + J2 


= In 


Vho 

W 


+ Fi+ a 2 In 


W'o 

W 


+ F 2 


(C.70) 




where a = a{UV,UW) and F = F{UV,UW, fii) are functions. Finally we write Z in the original 
coordinates as follows 


■ A 

(. w \ 1 , / 


1 -) + -aln ( 

'yw 

V Wout J 1 ^ 

r A 


— 

^1 + aw Inliw) + wGj 

'yw 


\ w J 


F 


Z = Z exp(/o + Ii) 
= Z exp 

= Z exp 


where a = a{uv \^and G = G{uvl^^' , fXi) are G^ functions with fXi = ViW 


(C.71) 


The transition 


In this case the time of integration is given by T = In ( ^ ) . Such a time of integration is obtained from 

the equation v', = P^Vj. The we have 


( 


i/ft 


u = u \ — 


Vj 


= w(^\ 

^ \VjJ 


1//3 


ft/ft- 


7/ft 


(C.72) 


w 


It then only rests to compute Z. Following similar arguments as for the transition II™ we get in this case 

A 

[ jw 


Z = Z exp 


^I + a'w ln(r'j) + w& 


(C.73) 


where now 

a' = a'{uvy^\uw^^'^) 

G' = & [uv]^ , Hw, h^i) 

are G^ functions with = wv^ and /ii = ViVj 

Saddle 2 

In this case we investigate the transitions of a vector field of the form 

{ u' = wu 

V'. = -PjWVj 

/ 2 

W = —^w 

Z' =-gZ, 

48 


(C.74) 

□ 


(C.75) 



where the coefficients Pj, 7 are positive. We assume that u G M'*'. Observe that now, in contrast with case 
1, we only have one expanding direction, which is u. This makes the study of the transition easier. Due to 
the same reason, it is more convenient to study a transition 


n“ : Sen Sex, 


(C.76) 


where to be general, we let Sen be any codimension 1 subset of obtained by setting one of the 

coordinates {v,w) to a constant and with u < Uouu and where 


Sex — 


u,v,w,Z)\u = Uout} • 


(C.77) 


Proposition C.5. Consider the vector field Y given by (C.75) and let Sen, Sex and 11“ be as above. Then 

0. 


f u \ 

v^= Vi [ - 

\ 'Uout / 

~ ( “ V 

w = w \ - 

\ Uout ) 

Z = Zexp — —— 1 + arc ln(w) + rcP 
[ 7 W \\ u ) 

where a = a{u^^Vi,u0w) and F = F{u^'Vi,u'^w,u) are functions. 

Proof of proposition C.5. We have that the time of integration is T = In 

f u Y' 

Vi=Vi[ - 

\UontJ 

f n Y 

w = w I - 

\UontJ 

It only remains to compute Z. Following similar arguments as in case 1 we have 


Z = Z exp 



{ llout \ 

jw \ 

. V u / 


1 + aw ln(t 6 ) + wF 


where a = a(u^’Vi,u'^w) and F = F{u^^Vi,u'^w,u) are functions. 


(C.78) 


It follows that 


(C.79) 


(C.80) 

□ 
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